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We study the Drell-Yan process mediated by a new bosonic resonance at the LHC. The bosons 
of spin-0, 1, and 2 with the most general leading-order couplings to Standard Model fermions and 
giuons are considered, which provide a model-independent formulation for future exploration of 
the resonance properties, such as its spin, mass and couplings. In the case of neutral resonances, 
we demonstrate how the shapes of the kinematical distributions change as one varies the chiral 
couplings of the quarks and leptons, and show how to analyze the couplings by making use of the 
forward-backward asymmetry. In the case of charged resonances, we propose a novel technique to 
effectively reconstruct the angular distribution in the center-of-mass frame, to a large extent avoiding 
the two-fold ambiguity due to the missing neutrino. Similar to the case of a neutral resonance, the 
spin information of the resonance can be extracted unambiguously, and chiral couplings and the 
asymmetries can be explored in a statistical manner. With the current LHC data, we present 
bounds on the mass and cross section times branching fraction of the new resonance and estimate 
the future reach. 



I. INTRODUCTION 

With the start of the LHC experiments, we have en- 
tered a new era of high-energy physics that directly 
probes Nature at the TeV scale. Depending on the un- 
derlying theory, new particles of different kinds may lead 
to novel and distinctive signatures at the collider. In 
addition to the highly anticipated discovery of the ori- 
gin of electroweak symmetry breaking, we will likely also 
discover other new resonances associated with this scale, 
presumably through the classic Drell-Yan (DY) process 
with a striking signal at around the TeV scale. If such 
a new resonance is indeed observed, it is important to 
determine many of its properties, such as the spin J, 
fermionic chiral couplings, and so on, in addition to its 
mass, width and electric charge. In this work, we con- 
sider new spin-0, 1, and 2 resonances that can contribute 
to the s-channel Drell-Yan processes^. In view of many 
possible models predicting such resonances, we keep the 
couplings between the standard model (SM) fermions and 
the resonances as general as allowed by symmetries. 

For a new neutral resonance, the final state involves 
two charged leptons whose momentum information can 
be fully registered by the detector. Therefore, with a 
sufficient rate, the particle mass can be readily deter- 



^ There have been many studies of resonances that can contribute 
to the Drell-Yan process. See, for example, [ll 456ll , where [3 435^ 
deal with a neutral boson, [3h445| deal with a charged boson and 
[HJildeal with both. 



mined from the invariant mass distribution. By boosting 
to the centcr-of-mass (CM) frame of the two leptons and 
noting that the boost direction preferably coincides with 
the momentum of the colliding quark, one can study the 
angular distribution of the leptons to extract the spin 
information of the resonance. 

The case of a charged resonance is more complicated. 
This is because the final state contains a charged lep- 
ton and the associated neutrino, thus missing energy is 
involved in such events. The resonance mass can be 
best determined from the Jacobian peak in the transverse 
mass distribution. However, there is a difficulty in find- 
ing the correct CM frame of the charged lepton and the 
neutrino. Even if one assumes that the resonance mass 
has been measured, there are generally two possible solu- 
tions for the longitudinal momentum of the neutrino. We 
develop a novel technique to effectively construct the an- 
gular distribution for the charged lepton in the CM frame 
in a statistical manner, and show how the spin and chiral 
couplings of the resonance can be extracted in a similar 
fashion as in the case of neutral resonances. 

This paper is organized as follows. In Section [Hi 
we classify types of new resonances that can mediate 
s-channcl Drell-Yan processes at the LHC, and discuss 
some of the general features of these processes. In addi- 
tion, we discuss the most general interactions between 
SM fermions and the new resonances allowed by the 
Lorentz and gauge invariances, with the corresponding 
Feynman rules given in Appendix \^ Section IIIII is de- 
voted to the discussions of the Drell-Yan process me- 
diated by a neutral resonance. The Drell-Yan process 
associated with a charged resonance is analyzed in Sec- 



tion llVI The current results of the LHC are used to place 
a bound on the masses and couplings of these new reso- 
nances in Section |Vl The findings of this work are sum- 
marized in Section I VII Appendix |B] contains details of 
the FeynRules implementation of these new interactions. 
Appendix C gives a short review of the Wigner d^^ „j, 
functions that are useful for the helicity amplitudes of 
our calculations. 



II. GENERAL INTERACTIONS OF NEW 
DRELL-YAN RESONANCES 



We will concentrate exclusively on color-singlet 
neutral-current and charged-current s-channel reso- 
nances contributing to the Drell-Yan processes at the 
LHC 
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A. Spin-0 states 

We begin with the neutral scalar boson S which can 
have the following Lagrangian 



C,, 



fi {9L 



^9k,l^)f,S-\fF;^,F^'^-^S, (2) 

where / can be either a quark or a lepton. The indices i 
and j run over generations, and the generation matrices 

9sij ~ i9s)ij ^^'^ 9pij ~ i9p)ij ^'"^ required to be Hermi- 
tian by the Hermiticity of the Lagrangian. F" denotes 
the gluonic field strength tensor. Here and henceforth, 
A denotes the cutoff scale of the effective interactions, 
which should be at least at the order of the resonance 
mass or higher. For the charged scalar boson 5^, we 
have 



C 



s± 



where hl;^. = {h{)^j and hp^^ 
general complex matrices. 
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p)ij are allowed to be 



where i generically denotes either an electron or a muon, 
V denotes a neutrino or an antineutrino and X the inclu- 
sive hadronic remnants. We consider the most general 
couplings for leading-order operators allowed by certain 
symmetries for such a new particle. This particle must be 
a color singlet and have integer spin. We concentrate on 
the possibilities of a scalar (S), vector {V) and traceless 
symmetric second-rank tensor (T), although other spins 
are possible in principle. There are two cases for the elec- 
trical charge of the boson. In the case of neutral-current 
processes, this boson (generically denoted by Rq) must 
be neutral whereas for the charged-current processes, this 
boson (generically denoted by Re) must have charge ±1. 
In Table U we summarize the properties of these new 
bosons along with the processes to which they contribute. 



Notation \Qe 



J 



Partonic processes 



Ro 0,1,2 uu,dd,gg ^l+r 

Re 1 0,1,2 ud^£+iy,du-^rv 

TABLE I: Resonance particles, their quantum numbers, and 
s-channel Drell-Yan processes. Qe and J represent their elec- 
tric charge and spin respectively. 



We now write down the most general Lagrangian be- 
tween these new bosons and the SM fermions and glu- 
ons allowed by Lorentz, quantum chromodynamic and 
electromagnetic gauge invariance. In each case, we 
only include the leading effective terms which are ei- 
ther dimension-4 or dimension-5 operators. Furthermore, 
we drop all terms which vanish when the masses of the 
initial-state and final-state particles are taken to zero. 



B. Spin-1 states 



Using the same notation, the most general Lagrangian 



for the neutral vector boson y^ is 



^v = fn' 



9nj + 9mj^5) fjV^, 



(4) 



where g(,^- = {gy)^J and 5^-^- = {9^)13 are required to be 
Hermitian matrices.^ We have dropped interactions of 
the neutral vector boson V with gluons since they do not 
contribute to this Drell-Yan process. The charged vector 
boson has the Lagrangian 



Cyi - un'^ [hl,^ + /4.,75) d,v+ + h.C. 

+^^^r {hUj + h^Mfl^) i^y^ + h.C, (5) 

where hy^- = ihy)ij and /i^- ■ = (/i^)ij a-i'c allowed to be 
general complex matrices. 



C. Spin-2 states 

The neutral tensor Lagrangian is given by 



^- = A 



/, {9L3 - 5iT.,75) irO'^f, + Yd^f,) 



^ Our convention is fixed with respect to the SM coupUngs as 
9zy = 9z{\tI - Qt'sl,), g{j^ = -^T^gz, where gz = g/c^, 
g is the weak coupling and Su, and c^ are the sine and co- 
sine of the Weinberg angle. The left- and right-chiral cou- 
plings arc related to the vector and axial-vector couplings as 
9RxL = QRxV — gUxA, gRxR = 9RxV + gRxA- 



- (drf^i- + d^f.^n (.41- + ff^kTs) /. 



T, 



flU 



(6) 



where the coupHngs g^^^ = {g^)ij and g^^ij- = iaAr)^] 
are general 3x3 complex matrices. Note that a 
dimension-4 operator between SM fermions and neutral 
tensor particle is not allowed since it would be propor- 
tional to the trace of the tensor which we have assumed 
to be traceless. Finally, the charged tensor has the inter- 
action Lagrangian as follows 



Ct± ~ -r 



A 






^AT»j75 ) dj 



Tj 



d^V.r) {I 



Tij 



hATijls) £j 


T+ 


(7) 


^J^ ^Ty = 





+h.c., 

where /i^.^. = (/i^)y, h{p^j = ih{^)i 

and hy^xii ~ (^yiT)y ^^^ general complex matrices. The 
corresponding Feynman rules have been worked out and 
can be found in Appendix [Xj 

III. NEUTRAL BOSON RESONANCES 

Diagrams for the Drell-Yan process pp — >■ £'^£~X medi- 
ated by a neutral boson at tree level are shown in Fig. [TJ 
including the SM diagrams with j/Z exchanges. The 
corresponding helicity amplitudes are listed in Table HIl 
where 



D 



X 



Ml 



iMxT 



x^ X- 



(8) 



We have expressed the scattering amplitudes in terms of 
the Wigner d^ ^, functions, where j is the total angu- 
lar momentum and m and m' are the difference of the 
hclicities of the initial-state and final-state particles, re- 
spectively [531 ■ ^ short review of the d^ ^, functions can 
be found in Appendix[Cj For s-channel scalar resonances, 
the initial state particles must have the same helicities to 
conserve angular momentum. The same is true for the 
final state particles. For this reason, we find that only 
the helicity combinations (A, A) — > (A', A') are nonzero. 
On the other hand, for s-channel vector and tensor inter- 
actions, angular momentum is not sufficient to determine 
the helicity combinations of the external particles. How- 
ever, the combined effect of the masslessness of the ex- 
ternal fields and the properties of the interaction vertices 
only allows opposite helicities for the incoming particles 
and also for the outgoing particles (A, —A) — >■ (A', —A'). 



A. Invariant mass spectrum 

The best observable to discover a new neutral boson 
coupling to quarks and charged leptons is in the spectrum 
of the invariant mass 



Mu = {Pe+ +Pe-) 



(9) 



If there is no interference, the shape of the invariant mass 
distribution is of a Breit- Wigner form and peaked at the 
mass of the new boson. However, if there is significant 
interference between the new resonance and the SM di- 
agrams, there can be appreciable changes in the shape. 
In particular, the peak may even be shifted. We show 
the di-lepton invariant mass distribution in Fig. [2] for 
the process pp — >■ £'^£~ X , only including the subprocess 
uu -^ £^£~ , mediated by a scalar (green dot-dashed), 
vector (red dotted) and tensor (blue dashed) boson, along 
with the SM (black solid) contributions including the full 
spin correlations. To see the interference effects clearly, 
we have adjusted the resonance rate to be the same value 
as the SM background rate near the peak. The row and 
column headers specify the nature of the chiral couplings. 
(The scalar fields have a factor of i accompanying the 75.) 
Here and henceforth, for illustration, the mass of the new 
particle is taken to be 1 TeV while the width is taken to 
be 20 GcV. We adogt the parton distribution functions 
(PDF's) CTEQ6L ^. The LHC energy is set at 7 TeV 
unless stated otherwise. 

For massless fermions, the scalar particle amplitude 
{M-s) is nonzero only for initial and final states of the 
same helicity while the SM contribution is nonzero only 
for initial and final states of opposite helicities. For this 
reason, there is never any interference between the two, 
and a scalar field always renders a Breit- Wigner shape 
peaked at the mass of the scalar particle. This can be 
seen analytically in Table |TT1 and also numerically in the 
dot-dashed (green) curves of Fig. [2] 

The new vector boson does interfere with the SM. The 
amount of interference depends on the parity properties 
of the couplings. The interference always flips sign at the 
mass of the new boson due to the phase change in the 
propagator s — My. In addition, there is an overall sign 
coming from the couplings. After summing over helicities 
and integrating over the scattering angle 0, we find that 
the interference with the SM photon diagram is given by 

^ / dco^e{M^M*y + M;Mv) (X 7^ {g'y„,,g'vjn) , 

hel •'^^ 

(10) 
where TZ{x) means the real part of x. The interference 
with the photon is independent of parity violation and 
only depends on the sign of the vectorial coupling. The 
interference with the Z boson diagram, on the other 
hand, is given by 

Y, f dcose{MzMlr + M*zMv)(x (11) 

hel 

'^ [[9zv9v7m + 9zA9Ami) \9zv9vjn + 9zA9Ajn)\ J 




FIG. 1: The s-channel Feynman diagrams responsible for pp -^ t^l^ X. Contributions to qq -^ t^l^ are from a photon (7), a 
Z boson {Z), a new scalar particle (S), a new vector particle {V) and a new tensor particle (T). Contributions to gg — > t^i~ 
are from a new scalar field [gS) and a new tensor field (gT). 



which has more complicated dependence on the vector 
and axial vector couplings. Since g^y is very small, the 
sign of the interference is more strongly dependent on 
the axial coupling to leptons. Again, the interference 
can be seen analytically in Table |ll] and numerically in 
the dotted (red) curves of Fig. [31 

In the case of the tensor field, the term proportional to 
'^m m' does not contribute to interference in the invariant 
mass distribution since d^ ^, is orthogonal to d]^ „^, and 
this interference vanishes after integration over 0. The 
'^m.rn' term of the tensor amplitude is due to the off-shell 
effects and does not contribute to the peak at Mu ~ 
My- Consequently, the interference is very weak and 
not readily observable in the continuum invariant mass 
distribution. The final result is that the tensor boson 
has a Breit-Wigner shape peaked at its mass as can be 
seen in the dashed (blue) curves of Fig. [2l which is hardly 
distinguishable from the case of a scalar. 



B. Transverse momentum distribution 

Once a new boson resonance is established in the in- 
variant mass spectrum, it will be of ultimate importance 
to study its other quantum numbers, such as its spin, chi- 
ral couplings, etc. For a massless particle, the transverse 
momentum is related to the scattering angle 



Pt = E sin ( 



(12) 



Thus, the differential distribution of px may contain ad- 
ditional information of chiral interactions of the reso- 
nance via the interference with the SM diagrams. We 
present the transverse momentum distribution of the neg- 
atively charged lepton for a variety of parity violation 
cases in Fig. [31 As expected, there is no interference 
between a spin-0 state and the SM diagrams. The fact 



that a spin-2 resonance interferes with the SM in the 
transverse momentum distribution gives us a potential 
new way of determining the spin of the resonance, unlike 
the invariant mass distribution. If we find interference 
present in both the invariant mass distribution and the 
transverse momentum distribution, we can conclude that 
the new resonance is a spin-1 particle. Furthermore, the 
amount of interference and the sign of the interference 
can give information about the size and sign of the parity 
violation in the couplings. We emphasize that the anal- 
ysis of the pt distribution does not require a knowledge 
of the quark moving direction, nor the reconstruction of 
the CM frame. 



C. Angular distribution 

One of the main advantages of the present process is 
the feasibility to fully reconstruct the CM system of the 
two charged leptons that is the rest frame of the new bo- 
son. Although we do not know the direction of the quark 
on an event-by-event basis, it is strongly correlated with 
the direction of the CM frame of the charged lepton pair 
due to the parton distribution functions of the quark ver- 
sus the anti-quark in a proton [46|. We calculate this 
angle by first boosting into the CM frame of the charged 
leptons, and then taking the angle between the moving 
direction of the negatively charged lepton and the direc- 
tion of the boost. In the case of the gluon initial state, 
both directions are equally valid and we simply use the 
direction of the boost to measure the negatively charged 
lepton as in the quark case. 

The angular dependence comes from the Wigner d^^ ^, 
functions. For each helicity combination and each spin 
for the new boson, these are determined by the kinemat- 
ics (see Appendix[C]). The mixture of the d^^ ,„, functions 
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TABLE II: Helicity scattering amplitudes for the parton-level processes. Tiie amplitudes correspond to the diagrams in Fig. [T] 
The particles in the s-channel exchange are labelled by subscripts (7, Z, S, V and T) while gS and gT indicate that the initial 
states are gluons with scalar and tensor exchange, respectively. A and A' are the helicities and we define 'H± (M) — M ± M^ . 



encodes the information of its spin and chiral interac- 
tions. For the scalar field, only cIq g contributes and so 
the angular distribution is flat. This can be seen in the 
dot-dashed (green) curves of Fig. |4l For a vector field, 
the rf±i ±1 functions contribute. Each one is squared and 
summed with the appropriate factors (see Table HI)) . This 
gives the angular distribution 



that of the axial coupling, as in 1 ± 75). However, it is 
interesting to note that the distributions are identical for 
either purely right-handed or purely left-handed, and are 
parity-transformed for right- and left-mixed couplings. 

For the tensor field, the angular dependence is mainly 
a mixture of the d'j_i ^.^ functions (the d^.^ j_i functions 
give a small perturbation) 
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Av{l + cos ef + Bv{l~ cos ef 



{9lJ.S + {9l,9Lf- 



(13) 



This is the formula of a parabola versus cos where the 
amount of parity violation determines where the mini- 
mum lies. This parabolic shape can be seen in the dotted 
(red) curves of Fig. |4l We find that the angular distribu- 
tion is symmetric whenever either coupling is pure vector 
or axial vector. If both the quark and the lepton cou- 
plings are parity- violating (mixtures of 1 and 7,5 terms), 
then the angular distribution is shifted to one side or 
the other. The shift is maximal when both couplings are 
purely chiral (the magnitude of the vector coupling equals 



Y,\Mj 
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Bt {I ~ cos ef {2 cose + if 



where At and Bt are of the same form as in Eq. (|13p . 
but with the tensor couplings. This is a quartie function 
with one local maximum and two local minima, a "W" 
shape. It can be seen in the dashed (blue) curves of Fig.|4l 
In particular, we find that if either of the couplings are 
purely tensor (g ex 1) or purely axial tensor {g ex 75) 
the distribution is symmetric, with the local maximum 
occurring at cos = 0. However, if both quark and lep- 
ton couplings are mixtures of 1 and 75, then the curves 
shift toward one side or the other, leading to parity vio- 
lation. As in the vector case, the shift is maximal when 
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FIG. 2: (Color online) Invariant mass distribution for the process pp —t- l^i~ X. The vertical axis is the differential cross 
section in arbitrary units and the horizontal axis is the dilepton invariant mass running from 950 GeV to 1050 GeV. The row 
and column headers specify the nature of the chiral couplings. The solid (black) curve is for the SM, the dot-dashed (green) 
curve includes the scalar field, the dotted (red) curve includes the vector field, and the dashed (blue) curve includes the tensor 
field. The scalar and tensor curves are indistinguishable and right on top of each other. The mass of the new particle is taken 
to be 1 TeV while the width is taken to be 20 GeV. The LHC energy is set at 7 TeV, and the CTEQ6L PDF sets are used. 
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the coefficient of 1 has the same magnitude as that of 75 . 
To be more reahstic, we fold in the CTEQ6L parton 
distribution functions and adopt some basic acceptance 
cuts on the transverse momentum and its pseudo-rapidity 
for the charged Icptons 



PTi > 20 GcV, \r]i\ < 2.5. 



(14) 



We see the effects in the three panels on the right-hand 
side in Fig. U) Due to the misidcntification of the correct 
quark momentum direction, the far forward/backward 
regions are diluted. The rapidity cut limits the angle 
reach in the same region as seen at the drop near 0.8. 



Co = will lead to the largest event rate, while a particu- 
lar choice of Cg may optimize the size of the asymmetry^ . 

It is clear that a scalar field will give ApB = for any 
value of cq since it has a flat angular distribution. 

For a vector field, it is natural to use Cq = 0. It can 
be shown that the forward-backward asymmetry at the 
peak of a vector resonance is given by 
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(16) 



[{9v) + (ffA) ] [i9v) + (5a) J 4 



D. Forward-backw^ard asymmetry 



where (jjr defines the ratio of vector and axial couplings 



It is customary to construct the forward-backward 
asymmetry, based on the partially integrated rates in the 
opposite angular regions, which can be defined as 



ApBico) 



N{coa0 > Co) - N{cos0 < -cq) 
N{coae > Co) + N{cose < -cq) 



(15) 



^ The asymmetry can be defined with any differential form with 
respect to the angular range. Other Lorentz invariant asymme- 
tries have also been proposed, for example 1591. I6OII . Others have 
noted that the asymmetry can be used to discriminate between 
different Z' models [3I. [sl. CT [sH [elMel . 
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FIG. 3: (Color online) Transverse momentum distribution of ^^ for the process pp — >■ ^+^^X. The vertical axis is the differential 
cross section in arbitrary units and the horizontal axis is the transverse momentum running from 450 GeV to 525 GeV. The 
row and column headers specify the nature of the chiral couplings. The curve legends are the same as in Fig. [2] 



and is given by 
cos (j)f - 
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For example, (f>f = corresponds to a pure vector cou- 
pling with g^ (X 1, (pf = 7r/4 corresponds to a pure chiral 



coupling with g' ex 1+75, and so on. Eq. (|16p shows 
that the absolute value of the asymmetry is bounded to 
be less than or equal to 3/4. Also, if any of the cou- 
plings gy, g\, gy and g^ is identically zero, there is no 
forward-backward asymmetry, resulting in parity conser- 
vation. The SM asymmetry for the uu — > ^+£~ process 
at ^/s 3> Mz is approximately 0.6. As mentioned be- 
fore, the quark momentum direction is ambiguous in pp 
collisions. We thus take the angle in the CM frame with 
respect to the boost direction, which is more likely to be 
the direction of the quark. 

Again, to be more realistic, we convolute with the par- 
ton distribution functions, and take the acceptance cuts 
as in Eq. p^ . In Fig.[SJa) and[5ljb), we present the equal- 
valued contours of forward-backward asymmetry for a 
vector resonance in the 4>X.-4>i plane. Fig. O^a) shows 
contours for uu — > l~^(.~ via the resonance only, and Fig. 
[5]Jb) includes the PDF convolution at the hadronic level. 



We see that the asymmetry is reduced when the parton 
distribution functions are turned on and acceptance cuts 
imposed. The SM value becomes ^ 0.37. As for the sig- 
nal, the maximum is reduced from 3/4 to ^ 0.41. The 
reason for this is that the boost direction corresponds 
with the quark direction much of the time but is some- 
times in the opposite direction. The asymmetry breaks 
up into quadrants separated by the angles (f)^ = rm /2 for 
some integer n where the asymmetry is 0. This is where 
the sign of the vectorial or axial coupling changes sign, 
thus changing the sign of the asymmetry. 

For a tensor particle, we find that cq 7^ is required 
to obtain a non-zero asymmetry. This can be seen by 
noting that for co = 0, 
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for all values of At and Bt- In other words, although 
the tensor distribution is asymmetric, it always shifts in 
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FIG. 4: (Color online) Angular distributions of £~ in the CM frame for the process pp — >■ £^£~ X. The curve legends are the 
same as in Fig. (2] In the top row purely scalar, vectorial and tensorial couplings are used. In the second row left chiral couplings 
are used while in the third row left chiral couplings are used for the quark but right chiral couplings are used for the leptons. 
In the left column the PDF's are turned off and the lab frame is the same as the CM frame, while on the right, the PDF's are 
turned on and the boost direction is used. On the right panels, the visible bend-down feature in the forward-backward regions 
is due to the kinematical cuts as described in the text. 



such a way as to have equal area under the angular dis- 
tribution for positive and negative values of cos 6*. For 
this reason, it is necessary to use cg 7^ 0. We find that 
the difference is maximized when cq = l/\/2- 
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(17) 



where (/)T defines the ratio of tensor and axial ten- 
sor coupHngs, similar to the vector case. Similar to 
Figs.[5lja) and[5jb), the equal- valued contours of forward- 
backward asymmetry for a tensor resonance are pre- 
sented in FigsjSJc) andlSfd). We again see that the asym- 
metry is reduced when the parton distribution functions 



are turned on. The maximum asymmetry is reduced to 
a little greater than ^ 0.37. The sign is similar to the 
vector case. 



If an experiment discovers a resonance in the invariant 
mass distribution of the neutral DrcU-Yan process, mea- 
suring the forward-backward asymmetry could help dis- 
tinguish the spins in the early days before enough events 
are accumulated for reconstructing the angular distribu- 
tion. If a significant asymmetry is found in Aps (0) , this 
would be evidence of a new vector resonance. If no asym- 
metry is found in Apb{0), then one should further deter- 
mine Afb{^/V2). li ApBi^/V^) is significantly different 
from zero while Afb{0) = 0, then this would be evidence 
of a new tensor particle. If both are zero, then this is evi- 
dence for a symmetric distribution, but could come from 
any of the spins we have discussed here. 
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FIG. 5: Contours of forward-backward asymmetry labelled by the values of Afb in the uu -^ t^l~ process. In the first row 
(a) and (b), a vector resonance is presented and in the second row (c) and (d), a tensor resonance is presented. Plots (a) 
and (c) show the result with the new physics amplitude alone in the lab frame while plots (b) and (d) include the PDF's and 
acceptance cuts. 



IV. CHARGED BOSON RESONANCE 



Transverse mass distribution 



Charged bosons contribute to the Drell-Yan process 
pp — >■ (.^vX through the diagrams in Fig. |6]at tree level. 
The resulting amplitudes are shown in Tabic IIIIl 



The first place we will look for a charged boson in 
Drell-Yan processes is in the transverse mass. Since we 
can not detect the neutrino, we can not fully reconstruct 
the invariant mass of the lepton system. The best one 
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FIG. 6: The s-channel Feynman diagrams responsible for pp -^ £vX. Contributions to this process are from a W boson {W), 
a new scalar particle (5), a new vector particle (V) and a new tensor particle (T). 
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TABLE III: Helicity scattering amplitudes for the parton-level processes. The amplitudes correspond with the diagrams in 
Fig. [6] The particles in the s-channel exchange are labeled by subscripts {W , V, T and S). X and A' are the helicities and we 
define T±{M) = M ± M where M means that we replace /i by /i in M. 



can do is to construct the transverse mass which contains 
the charged lepton momentum and the missing transverse 
momentum: 



Mj^ = {EtI + ETmiss) — {VTt + PTmiss) 



(18) 



where the transverse energy is defined as Et = 
\J M'^ +Pt' ^^'^ PTmiss is identified as pTi/- In practice, 
we assume that the SM Icptons are massless. Further- 
more, in a 2 — > 2 process in the absence of transverse 
motion, the missing transverse momentum is equal and 
opposite to that of the charged lepton. This allows us 
to simplify the transverse mass to Mt = 2pTe- This 
transverse mass distribution develops a Jacobian peak at 
the mass of the resonance particle. We plot examples 
of transverse mass distributions for a variety of parity 
violation cases in Fig. [T] Once again, we see that the 
shape of the transverse mass distribution depends on the 
parity properties of the fermionic couplings^. Both the 



vector and the tensor shapes are modified, due to the 
interference effects near the resonance peak, while the 
scalar shape remains fixed. This could be important in 
determining the mass of a new charged resonance. More- 
over, the discovery of interference in the transverse mass 
distribution implies that the spin of the new resonant 
particle is greater than while the sign and size of the 
interference can give information about the sign and size 
of the parity violation in the couplings. We emphasize 
that the analysis of the transverse mass distribution does 
not require a knowledge of the quark moving direction, 
nor the reconstruction of the CM frame. 



B. Angular distribution 



* The importance of the transverse mass in determining the prop- 
erties of the couplings has also been discussed in |67|| . 



Unlike the case mediated by a neutral boson, the angu- 
lar analysis is known to be difficult for the charged boson 
mediation due to the missing neutrino in the final state, 
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FIG. 7: (Color online) Transverse mass distribution for the process pp — >■ H vX. The vertical axis is the differential cross 
section in arbitrary units and the horizontal axis is the transverse mass running from 900 GeV to 1050 GeV. The row and 
column headers specify the nature of the chiral couplings. The curve legends, mass, width, energy and pdf set are the same as 
in Fig. [3 



especially for the LHC as a symmetric pp collider^. As- 
suming that pti/ = — Pt£, we can solve for pzv in terms 
of the measured charged lepton momentum and the lep- 
tonic invariant mass Mi^, 



Vzl 



( Mf 



Iv 



\2p\ 



re 



-1 ± 



pre 



IMI 



(19) 



where all quantities are in the lab frame. If the width 
of the new resonance is sufficiently small, then the in- 
variant mass is well approximated by the resonant mass 
Mi^ ~ M. With a clear signal identification, the res- 
onance mass could be measured in the transverse mass 
distribution, just like the Mw determination in the SM. 
Even so, we still have a two- fold ambiguity in pz^. On 
an event-by-event basis, we do not know which one is 
correct. Instead, for each solution, we can calculate the 
angle of the charged lepton in the CM frame with respect 
to the boost direction to approximate the quark moving 



^ Other attempts to reconstruct the spin with a missing particle 
can be found in 168 



direction, formally defined by 



cos 



pCM 



(20) 



which in turn also suffers from the above ambiguity. We 
denote the smaller solution of \pz}y\ by the subscript S 
and the larger solution by the subscript L. The cosines 
of these angles in the CM frame, can be expressed by the 
lab quantities as 



cos 6*3 
cosOl 



4p; 



AP 



Ti sign [-—-E, 



M 
2 



1 



Mi 

M2 ■ 



(21) 



When El = M/2, this corresponds to the situation where 
the lab frame and the CM frame coincide. Both solutions 
of cos6' are identical when Ei < M/2 ^. On the other 



® However, there is still the possibility that the boost direction does 
not coincide with the quark direction as in the neutral current 
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hand, the two solutions differ by a sign when Ei > M/2. 
In this latter case, when a wrong solution is chosen, it 
simply moves the angle from cos 9 — > — cos0. We show 
a contour plot of cos 8 in the p^i — pre plane for the 
two solutions in Fig. |8l As is clear in the plot, cos 9 
does not depend on p^e except at the transition point 
Ei ~ M/2 — \/p\i + p^i , where it flips sign. 
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FIG. 8: Contour plot of cos 6 as a function of the momentum 
of the charged lepton for the process pp — >■ lu. The horizon- 
tal lines are the contours. The contours are (starting from 
the bottom) ±0.99, ±0.90, ±0.75, ±0.50, ±0.25 and 0. The 
large solution is negative everywhere while the small solution 
changes sign at the transition point Eg — Af/2 indicated by 
the circular curve. 

These observations lead us to reconstruct adequate an- 
gular variables to compensate the loss of information due 
to the missing neutrino. By taking the large solution 
cos 9]^ , which is negatively definite, we obtain the distri- 
bution — I cos9\. Convoluting with the PDF and impos- 
ing the acceptance cuts 

pre > 250 GcV, pTrmss > 250 GcV, \r]e\ < 2.5, 

. (22) 
we reconstruct the angular distributions as seen in the 

left column of Fig. |9l We further symmetrize the dis- 
tribution by splitting the large solution in half and tak- 
ing the mirror image on the cos 9 > side, as seen in 
the right column of Fig. [9] with the solid (black) curves. 
With this prescription, we reproduce the correct angular 
distribution if the true distribution is symmetric {e.g., 
no parity violation), and we only obtain the average over 
the positive and negative regions of cos6', as compared 
with the dashed (blue) curves in the right column that 
present left-left chiral couplings for a vector and a tensor 
state. In all cases, the spin information is well preserved 
and this gives us an unambiguous way to determine the 
spin of a charged boson resonance without the need to 
reconstruct the CM frame. 

On the other hand, the small solution carries both signs 
and thus contains information not only about the spin, 
but also the asymmetry. In the region Eg < M/2, this 
solution is the same as the large solution and thus the 
angle is typically correctly reconstructed. In the region 
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FIG. 9; (Color online) Angular distribution of £ for the pro- 
cess pp -^ Iv when only the large solution of the neutrino p, 
is taken. The cut-off feature in the forward-backward regions 
is due to the kinematical cuts as described in the text. 



El > M/2, however, either this solution or the one with 
an opposite sign could be right and we cannot determine 
it on an event-by-event basis. In each plot in the left col- 
umn of Fig. 1101 we plot the angular distribution using the 
small solution. The distribution is different for each spin 
and asymmetry and it generally covers the entire range of 
cos 6*. As expected, the small solution is essentially cor- 
rect about cos 6* < —0.6, but there is a clear deficit when 
cos 9 approaches 0~. The missing events are incorrectly 
assigned to the bins in — cos 6* > 0, which is also seen as 
an excess near cos — > 0+ . 

Fortunately, we are able to simulate the expectation 
and thus restore the distribution on a statistical basis. 
In each plot on the right column of Fig. [TUl the dashed 
(blue) curve is the expected angular distribution if we 
know the full momentum information of the neutrino. 
For cos6' > 0, we take the difference between the solid 
(black) curve on the left column and the dashed (blue) 
curve on the right column, to obtain the excess. We 
would like to move those events back to the bins in 
— cos 6* < to restore the original distribution. It is 
very important to realize that the fractional excess (the 
above excess divided by the solid (black) curve on the 
left column) turns out to be numerically the same for all 
of the scalar, vector and tensor resonances. Denoted by 
/ as a function of cos9, we obtain this simulated frac- 
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tional excess as a universal function for all spins, shown 
in Fig. 111! This observation leads to a powerful procedure 
for the restoration of the angular distribution: Given a 
data set, presumably like the solid (black) curve on the 
left column, we apply the fractional excess function to 
the cos^ > region bin by bin, then subtract this re- 
sult out from the data, and finally move the result to the 
region of cos 9 < for the correction. 

Plotted on the right column in solid (black) are the 
corrected angular distributions. In the case of symmetric 
distributions, the reconstruction works perfectly and the 
solid (black) curve completely coincides with the expec- 
tation of the dashed (blue) line. In the case of asym- 
metric distributions, the reconstruction is perfect on the 
edges but slightly off in the middle where the dashed 
(blue) curve can be seen. This is due to the more likely 
mismatch between the directions for the true quark mo- 
mentum and the boost. Nevertheless, it is sufficient to 
present features of the asymmetry. We reiterate that 
our correction for each distribution is independent of the 
resonance spin, although it depends on the parton distri- 
bution function as well as the mass of the new resonant 
particle. 



C. Forward-backward asymmetry 

The angular distributions for the large solution average 
out the events and thus lead to no asymmetry. The cor- 
rected angular distributions for the small solution as in 
the right column of Fig. [10] preserve the asymmetry prop- 
erty to a large extent. We could use these angular distri- 
butions in the same way as in the neutral current case, or 
we could use the small solution directly only in the region 
I cos^l -^0.5 without performing the cos 9 -(-> — cos 9 cor- 
rection. Thus similar analyses to Fig. [5] in the the neutral 
currents can be performed. This is a significant progress 
for the charged boson signal at the LHC as a symmetric 
collider. 
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V. SEARCHES FOR DRELL-YAN TYPE 
SIGNALS AT THE LHC 

A. Neutral current channel 

The CMS and ATLAS CoUaborations have studied 
dilepton events at 7 TeV in search of a heavy neutral 
boson [8J, |85[ . They present their results in terms of the 
production cross section times branching fraction as a 
function of mass, which takes into account the couplings, 
mass and width of the new particle for general spin. In 
this note, we extend this by projecting the expected 95% 
confidence level bounds to 10 fb~^ at 7 TeV and 10 fb~^ 
and 300 fb~^ at 14 TeV. We also place bounds on the 
product of couplings as a function of mass for each spin 
while using a standard width. 



FIG. 10: (Color online) Angular distribution of £ for the pro- 
cess pp -^ iu when only the small solution of the neutrino 
Pz is taken. Each row lists a case of distinct spin and parity 
violation. Each plot on the left gives the angular distribution 
directly obtained from the small solution. On the right, the 
dashed (blue) curve is the angular distribution with full mo- 
mentum information of the neutrino. The solid (black) curve 
is the angular distribution obtained from the corresponding 
distribution on the left by following the procedure explained 
in the main text. The cut-off feature in the forward-backward 
regions is due to the kinematical cuts as described in the text. 



We simulate the process pp — )• l^l~ X in the SM at 
tree level using the parton- level Monte Carlo package 
CalcHEP [8a, [83 and normalize our parton- level calcula- 
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FIG. 11: Fraction of events that should be removed from the 
cos 8 > side of the small-solution angular distribution and 
added to the cosO < side of the distribution, universal for 
a charged resonance of any spin. 



tions to the data by fitting to the two highest bins (from 
85 GeV to 95 GeV) from Fig. 2 of Ref. ^. Including the 



num.erical normalization factors 0.56 for the e^e~ chan- 
nel and 1.01 for the fJ-'^ fi^ channel, our results and the 
CMS data are given in Fig. [12] by the solid (black) curve 
and the (black) dots, respectively. Signals, including 
interference with the SM background, corresponding to 
new resonances of different spins, masses, and fcrmionic 
couplings were then done and superimposed on the plot. 
For simplicity and illustration purposes, we assumed that 
all the couplings were taken to be real in the numerical 
studies and that only the first two generations of fermions 
were taken into account. Results for more general cases 
are expected to be qualitatively the same. 

In numerical estimates of the signal events, the width 
of the new resonance was taken as 3% of its mass, in line 
with the assumption made in Ref. [SJ]. We considered 
the events in the invariant mass window of ±20% of its 
mass 
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-M < 
5 
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(23) 



where M is the mass of the new particle. The number 
of events for the e~^e~ and yU^/i^ channels were added 
{N = Nf. + Nf^) to give the total number of predicted 
events for the SM. This number of events was plugged 
into the one bin log likelihood [57J 



LL = 2 



TV In 



TV 



(24) 



where v is the number of events (for both the e~^e^ and 
//"•"/x" channels) expected in the SM plus the new boson. 
A value of LL = 4 was taken as the 95% confidence level 
and v was solved for. Taking into account the SM expec- 
tation and given an available integrated luminosity, this 



FIG. 12: (Color online) Invariant mass distribution of the 
charged leptons for the process pp — >■ f^'^ fj,~ X at the LHC 
with 7 TeV. The solid (black) curve is the SM expectation 
while the dot-dashed (green), dotted (red) and dashed (blue) 
curves are respectively for spin-0, -1, and -2 resonances with 
a mass of 750 GeV and a width of 22.5 GeV. The latest CMS 
data [8J] are superimposed in the figure as the (black) dots. 



is converted into a signal cross section times the branch- 
ing fraction for a given mass, and the results are plotted 
in Fig. US] for 7 TeV and 14 TeV. 

We wish to emphasize that our results here can be 
broadly applied to any resonant signal as outlined in the 
earlier sections. The production cross section is governed 
by the resonant coupling to the initial state partons, and 
the decay branching fraction is proportional to the cou- 
pling to the final state leptons. Thus with the determina- 
tion of the resonant mass and width by the kinematical 
peak and the line shape, we expect to gain the informa- 
tion for the fundamental couplings gqg£ for a scalar or a 
vector resonance, and gqg^/ h? for a tensor. 

To illustrate this point, we take the commonly studied 
"sequential Z' model" as an example for a vector reso- 
nance {Z'gj^j), which has the same coupling as the SM Z 
boson. In Fig. [131 '^e also include the cross section times 
branching fraction (the blue dashed curve) for this model 
as labeled by Z'gj^j. Where this curve crosses the solid 
(black) curve gives the bound (or projected bound) for 
this model. We find a lower bound on ^5^/ of 1135 GeV 
consistent with the CMS results. We also find that CMS 
could bound the mass at ~ 2.5 TeV with 10 fb^^ at 
7 TeV, - 4.1 TeV with 10 fb'^ at 14 TeV and ~ 5.5 TeV 
with 300 fb~^ at 14 TeV. This procedure can easily be ap- 
plied to other specific models. The couplings and widths 
simply need to be set and the mass scanned over to de- 
termine the cross section times branching fraction. 



15 



pp^/+/"X (7 TeV) 



10 



-1 _ 



& 10-2 

m 

X 



10" 



\ 




' 










1 1 . 


-^.,,_^^ \ 
















— — __^s___ 
















\ 

\ 


\ 


\ 










CMS 








\ 


\ 






: 


"^^ 




- 






\ 


•s 


\ 


\ 


-1 , , , , 1 , , , 














\ 



500 1000 1500 2000 2500 3000 
M (GeV) 



10 



-2 F 



m 

X 

b 



pp^rrX (14 TeV) 

' ' \ 




10- 



10" 



10- 



1000 2000 3000 4000 5000 
M (GeV) 

FIG. 13: (Color online) 95% C. L. bound on the production 
cross section times branching fraction versus the new resonant 
boson mass. The curve labeled CMS is with respect to the 
data measured by the CMS Collaboration [8J] , while the other 
solid (black) curves are projections for the specified integrated 
luminosity. The dashed (blue) curve is for a Z' with the same 
couplings as the SM Z. 



B. Charged current channel 

Following the same procedure as in the previous sec- 
tion, we calculate the 95% confidence level bounds and 
projected bounds on a new charged resonant boson. 
We base our results on the the CMS data published in 
[Ulsi]. We simulate pp -^ (.v {i = e, ^i) in the SM at tree 
level, again using the parton-level Monte Carlo package 
CalcHEP, and normalize our parton-lcvcl calculations to 
the data by multiplying our simulation with a numerical 
factor which brings the two highest bins (from 50 GeV to 
100 GeV) into agreement. We found the numerical fac- 
tor to be 0.76 for the electron and 0.91 for the muon. A 
plot of our SM calculation and the CMS data can be seen 
in Fig. [HI where the solid (black) line is the parton-level 



prediction of the SM and the (black) dots are the CMS 
data for the case of the muon. 

pp->//vX 
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FIG. 14: (Color online) Transverse mass distribution of the 
charged lepton at the LHC with 7 TeV. The solid (black) curve 
is the SM expectation while the dot-dashed (green), dotted 
(red) and dashed (blue) curves are respectively for spin-0, - 
1, and -2 resonances with a mass of 1 TeV and a width of 
20 GeV. The latest CMS data [8y| are superimposed as the 
(black) dots. 



The significance was calculated by taking the number 
of events in the window 



2 6 

-M < Mt < -M 
5 5 



(25) 



The one bin log likelihood was then calculated for the 
sum of the electron case and the muon case. 

We again plotted the LL = 4 line as the 95% confidence 
level in Fig. [15] Once again, to illustrate the approach to 
generalize to other models, we plot the cross section times 
branching fraction of a new charged vector boson with 
the same couplings as that in the SM, the "sequential W' 
model", denoted by Wgj^.^. We find that W'gj^j is bound 
to be heavier than 1500 GeV at 95% confidence level 
which is close to the CMS result [Sal. We also find that 
CMS could achieve a bound of ~ 2.5 TeV with 10 fb'^ at 
7 TeV, a bound of - 4.5 TeV with 10 ih'^ at 14 TeV and 
a bound of - 5.2 TeV with 300 ib'^ at 14 TcV. We again 
note that this procedure can be followed with any models. 
After the couplings and widths are set appropriately, the 
cross section times branching fraction can be calculated 
and plotted as a function of mass. 



VI. SUMMARY 

We have considered the most general new resonant s- 
channel contributions to Drell-Yan production of leptons 
at the LHC including spin-0, 1 and 2 bosons. Wc formu- 
lated the most general leading-order interactions between 
these new particles and the SM fields involved in the DY 
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FIG. 15: (Color online) 95% C. L. bound on the production 
cross section times branching fraction versus the new reso- 
nant boson mass. The curve labeled CMS is with respect to 
the data measured by the CMS Collaboration [83, [S^] , while 
the other solid (black) curves are projections for the specified 
integrated luminosity. The dashed (blue) curve is for a W' 
with the same couplings as the SM W . 



channel that satisfy Lorentz and EM invariances includ- 
ing both parity conserving and violating terms. Using 
these interactions, we have calculated the hclicity ampli- 
tudes and expressed them in terms of the Wigncr d?^ ^, 
functions, explicitly showing the angular dependence of 
these collisions in the CM frame. 

For the neutral current process pp — > £'^i^ X , we find 
that 

• the lepton pair invariant mass distribution may 
provide information for the chiral interactions. A 
new spin-1 field interferes with the SM process, 
thus modifying its shape from the Breit-Wigner 
resonance and the position of its peak, while both 
spin-0 and spin-2 fields do not have significant in- 
terference with the SM and appear as simple Breit- 



Wigner resonances. 

• the transverse momentum of the charged lepton 
also provides information for the chiral interactions. 
Both the spin-1 and spin-2 fields present interfer- 
ence with the SM background. Along with the in- 
variant mass distribution, this gives another way to 
distinguish between the spins. 

• defining an angle of £" with respect to the boost 
direction (likely to be the initial quark direction) 
in the CM frame of the system, one is able to con- 
struct the Wigner d^ ^, functions as well as their 
asymmetry due to parity violation. We note that 
an asymmetry only occurs for a spin-1 or 2 boson 
when both the quark coupling and the lepton cou- 
pling are not either purely vectorial (tensorial) or 
axial vector (axial tensor) and is maximized when 
the vectorial (tensorial) and axial (axial tensor) 
couplings are equal in magnitude. 

• although it is well-known that a chiral spin-1 bo- 
son generates an asymmetry when the forward and 
backward directions are defined in the full angular 
range from < cos 9 < 1 and — 1 < cos 6* < 0, re- 
spectively, the tensor asymmetry would be zero if 
taking this range. Instead, we find that the asym- 
metry is maximized for the tensor when considering 
the asymmetry in the range l/v2 < cos 9 < 1 and 
— 1 < cos 9 < — 1/-\/2, respectively. Thus the asym- 
metry gives a new way to distinguish between the 
spins. 

• scanning over mass and cross section times branch- 
ing fraction for the new bosons, we obtained the 
95% confidence level bounds based on the current 
CMS results as well as projected bounds for future 
integrated luminosities and machine energies. To 
illustrate how to apply our formulation and calcu- 
lation to other DY resonances, we plotted the cross 
section times branching fraction as a function of 
mass for a new vector boson with the same cou- 
plings as the SM Z boson (denoted by Z'gj^). 

For the charged current process pp — ?> i^vX^ we find 
that 

• the iv transverse mass distribution may provide in- 
formation for the chiral interactions. Spin-1 and 
spin-2 fields interfere with the SM process, thus 
modifying its shape from the conventional nonin- 
terfering transverse mass distribution as well as 
the position of the Jacobian peak while a spin-0 
field does not have significant interference with the 
SM and appears with the conventional shape. The 
shifted position of the peak can be important in 
determining the mass of the new resonance. 

• although determining the angle of i^ with respect 
to the boost direction is more challenging due to the 
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two-fold ambiguity in the z-componcnt of the neu- 
trino momentum, we find a novel statistical method 
for reconstructing the angular distribution. This 
method involves creating two distributions, one for 
the small neutrino z-component momentum and 
another for the large neutrino z-component mo- 
mentum. Wc find that the large solution faithfully 
preserves the symmetrized Wigner d^^ ^, function 
and therefore fully determines the spin of the new 
resonance but not the asymmetry. The small solu- 
tion distribution, on the other hand, does contain 
information about the asymmetry and we show how 
that information can be extracted. 

• scanning over mass and cross section times branch- 
ing fraction for the new bosons, we obtained the 
95% confidence level bounds based on the current 
CMS results as well as projected bounds for future 
integrated luminosities and machine energies. To 
illustrate how to apply our formulation and calcu- 
lation to other DY resonances, we plotted the cross 
section times branching fraction as a function of 
mass for a new vector boson with the same cou- 
plings as the SM W boson (denoted Wgj^j). 

We would like to reiterate that our formulation makes 
the future phenomenological and experimental searches 
straightforward. Our proposal for the large and small so- 
lutions for the charged current channel overcomes the dif- 
ficulty of the two-fold ambiguity due to the missing neu- 
trinos, that should be adopted for studies of new charged 
resonances. 



Appendix B: FeynRules Implementation 

We implemented our formulation using the FeynRules 
Mathematica package [90|- We will make the Feyn- 
Rules model files available through the FeynRules model 
database where those who are interested can download 
and use it. We briefly describe this implementation in 
this section. 

This implementation uses two FeynRules model files. 
The first is the SM that comes with FeynRules ("SM.fr"), 
although we will include a copy of this file in case changes 
are made to the SM model files later. We added all the 
new particles, parameters and Lagrangian terms to a new 
file we call "plus.fr" . 

The names of the particles along with their charge, 
mass and width are listed in Table IIVI These are the 





Name Mass Width 


s 


SV MSV WSV 


^M 


VV MVV WW 


J- IJ.U 


TV MTV WTV 


s± 


SVP± MSVP WSVP 


v^ 


VVP± MVVP WVVP 




TVP± MTVP WTVP 



TABLE IV: Particles implemented in FeynRules. In the first 
column is the symbol we use for the particle in this article. 
In the last 3 columns are the ASCII names we use for these 
particles, their masses and widths in FeynRules. These are 
the names that would be used in a simulation. The masses 
are set to 1 TeV by default while the widths are set to be 
20 GeV by default, but both of these parameters are free and 
can be set by the user. 
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Appendix A: Feynman Rules 

In this appendix, we provide the Feynman rules used in 
the analyses. The interaction vertices between the neu- 
tral resonances and the SM fields are listed in Fig. [Tni 
The interaction vertices between the charged resonances 
and the SM fermions are listed in Fig. [iTl The relevant 
SM vertices are reproduced in Fig. [TS] to fix our conven- 
tion. 



names that are used when running a Monte Carlo pro- 
gram. 

The couplings are implemented with names that are 
similar to the ones we use in this article. However, there 
is a six-character limit in CalcHEP for the names of pa- 
rameters, so some names are shortened. We list the pa- 
rameter names in Table [V] The couplings are imple- 
mented in a very general way in the FeynRules file and 
are all set to zero by default so that the user can turn 
on the vertices that they are interested in. Some are re- 
moved from the Monte Carlo file using "Definitions" in 
the FeynRules model file. These can easily be turned 
on by removing the appropriate definitions (the full cou- 
plings are included in the Lagrangians) , but the user 
should remember the hermiticity requirements on some 
of the couplings. Furthermore, if the user desires to use 
this implementation with a Monte Carlo package that can 
handle complex parameters, it is possible to turn off the 
splitting of these couplings into real and imaginary parts 
by commenting out the "Definitions" line of the complex 
couplings. 

The full Lagrangians described in this article were in- 
cluded in the FeynRules model files. FeynRules was then 
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s mh - iannriS"' 



S i(5s)/j - i9phl5 




s lUh - (gf^.n^S"' 



V, n"M)u + {9l)<j75]5'" 




d" 



V, n"[{9t-h + i9i)., 



75]5" 




V, ^Yligvhj + {gA).3^5. 







{glT)ijji]h^Pi+rp'i)S''' 
-ig%).MirP2+rp'2)S''' 



I'"' iKflf ),,■ - {g'AT\M (7"Pi + 7^pi') 

+il(ffr^),:.-(5^V)i,76](7''P2 + 7''P^) 






jf[-(9"°.9''" + S"'V")(Pi-P2) 

-{piP2+PiP^)g°'^+piP2g''^ 

+PlP2g'"^ + Pip2.9''" + PlP29'^]S° 





s 



s[(9r)u ■ 

+i[igPh ■ 



■{gAThi5]ii^Pi + Yp'i)S'' 
{g%hir,](YP2 + rp^2)S'' 



-Pi-P2g"''' + p"iP'i)S'" 



FIG. 16: Interaction vertices between the new bosons and SM fields. 



run on this file and all the resulting vertices were checked 
against our independent calculations of these vertices and 
agreement was found. 

The CalcHEP interface 19l| was then run on this 
model, which generated a set of CalcHEP model files. 
We used CalcHEP to generate the analytic formulas for 
the squared matrix element for all the 2 — )- 2 Drell-Yan 
processes where we only included the first generation 
fermions. We checked these formulas against our own in- 
dependent calculations of these squared matrix elements 



and found agreement. We then used these model files to 
do the numerical studies described in this article. 



These model files are intended to be used in unitary 
gauge and below the effective cutoff of these vertices. 
Furthermore, these interactions are not ultraviolet com- 
plete and should not be used for other processes than the 
Drell-Yan processes described here. 
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s+ [i(ft|),7 - {h^.n^^s^" 




s+ ^m„ - miji, 




4 



v; n'W^h + c^.Ms^'' 




+jVW)., - (h\r).M(iyi + YAW- 




v; n"\iH'h + (h'^hi,] 







FIG. 17: Interaction vertices between the new charged bosons and SM fermions. 




leQfY 



w+ f^yi^ - 75)% 




^/' *7''(flzv + 9zAli) 



w^;^7''(i-7om 



KM 



FIG. 18: Interaction vertices between the SM gauge bosons and fermions, where the weak couphng gw is related to the Fermi 
constant Gphy^ = ^. 



Appendix C: 






Function Review 



tides' hclicitics^. The final state has the same total angu- 
lar momentum and spin m' along the direction of motion 
of one of the final state particles and at an angle with 
respect to the z-direetion. We will define the x-direction 
such that the whole process occurs in the x-z plane. The 
spin along this direction is the difference of helicity of the 



In this appendix, we briefly review the Wigner d?^^ ^^, 
functions. Suppose the incoming state has total angu- 
lar momentum j and total spin m along the direction 
of motion of one of the particles. We will call this the 
z-direetion and the total spin is the difference of the par- 



"^ Although orbital angular momentum plays a role, it's z- 
component is 0. As a result, the total angular momentum z- 
component is equal to the spin 2-component (jz = m). 
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Name Real part Imaginary part 


akj 


gsfij 


gSfRij 


gSflij 


akj 


gPfiJ 


gPfRij 


gPfliJ 


9% 


gSg 






9U3 


gVfij 


gVfRij 


gvaij 


aLj 


gAfij 


gAfRij 


gAflij 


9% 


gVg 






9^,1^ 


gTfij 


gTfRij 


gTflij 


9'at.J^ 


gUfij 


gUfRij 


gUflij 


9'tI^ 


gTg 






^U 


hSfij 


hSfRij 


hSflij 


4. 


hPfiJ 


hPfRij 


hPflij 


<., 


hVfij 


hVfRij 


hVffij 


^A., 


hAfij 


hAfRij 


hAffij 


4., /A 


hTfij 


hTfRij 


hTflij 


4t,./A 


hUfij 


hUfRij 


hUflij 


/i^,,/A 


hYfij 


hYfRij 


hYffij 


^atJ^ 


hZfij 


hZfRij 


hZflij 



TABLE V: Couplings implemented in FeynRules. In the first 
column is the symbol we use for the coupling in this article. 
In the second column is the ASCII name we use for this cou- 
pling in FeynRules while the third and fourth columns contain 
the the real and imaginary parts. The letter "f ' refers to the 
flavor and runs over "u", "d" and "1" for the neutral cur- 
rent couplings while it runs over "q" and "1" for the charged 
current couplings. The i and j stand for the generation. 



final state particles. With these definitions, the d^^ ,^, 
function is defined as the overlap between the incoming 
and outgoing states 



C™'(^) = 0''"^'.^b''™> = (j'™' 



\e'-'\3.m) 



(Cl) 



where we have extracted the angular dependence into the 
operator exp{iJy0) which rotates the final state around 
the y-axis to an angle 9 with respect to the z-axis. Com- 
monly used functions are listed in Table I VII For j = 0, 



the generator of rotations around the y-axis is Jj, = 
and so 



For j = 1, the generator Jy is 



(C2) 




(C3) 



We also find that 




(C4) 



rfi — 

"0,0 — 


1 


di,±i = 


i(l±cos6i) 


di,±i = 


i(l±cos(9)(2cos6i=Fl) 


di.il = 


-i(l±cos6')sin6' 



TABLE VI: Wigner d-' , functions relevant for our processes 
and used in the helicity amplitudes given in Tables Ull and [Till 



and 



J, 



2)1 



Jf. 



j2n+l _ J 

y — J" 



(C5) 



where n > 0. With this information, exp{iJy9) can be 
expanded and it can be shown that 



/i(l + cos6i) ^sine* i(l-cose')\ 



dl 



JJy0 



V2 



V2 '"''' " """ " ^/2 "'" " 

\^i(l-cos6') ■^sin6' i (14- cose*) y 

. (C6) 
where m and m' refer to the elements of the matrix. 
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